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Abstract 

Event-by-event fluctuations of the net baryon number and electric charge in nucleus-nucleus 
collisions are studied in Pb+Pb at SPS energies within the HSD transport model. We reveal 
an important role of the fluctuations in the number of target nucleon participants. They strongly 
influence all measured fluctuations even in the samples of events with rather rigid centrality trigger. 
This fact can be used to check different scenarios of nucleus-nucleus collisions by measuring the 
multiplicity fluctuations as a function of collision centrality in fixed kinematical regions of the 
projectile and target hemispheres. The HSD results for the event-by-event fluctuations of electric 
charge in central Pb+Pb collisions at 20, 30, 40, 80 and 158 A GeV are in a good agreement 
with the NA49 experimental data and considerably larger than expected in a quark-gluon plasma. 
This demonstrate that the distortions of the initial fluctuations by the hadronization phase and, 
in particular, by the final resonance decays dominate the observable fluctuations. 
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I. INTRODUCTION 



The aim of the present paper is to study the fluctuations of the net baryon number and 
electric charge in nucleus-nucleus (A+A) collisions at SPS energies. We use the HSD [l| 
transport approach which reproduces both the different particle multiplicities and longitu- 
dinal differential rapidity distributions for central collisions of Au+Au (or Pb+Pb) from 
AGS to SPS energies rather well 2]. (see, e.g., Refs. , H El Q and 

references therein) reveal a new physical information. The fluctuations in A+A collisions 
are studied on an event-by-event basis: a given observable is measured in each event and 
the fluctuations are evaluated for a specially selected set of these events. We recall that the 
statistical model has been successfully used to describe the data on hadron multiplicities in 
relativistic A+A collisions (see, e.g., Ref. and a recent review |14,]) as well as in ele- 
mentary particle collisions This gives rise to the question whether the fluctuations, in 
particular the multiplicity fluctuations, do also follow the statistical hadron-resonance gas 
results. Recently the particle number fluctuations have been studied in different statistical 
ensembles Q]; the statistical fluctuations can be closely related to phase transitions in QCD 
matter, with speciflc signatures for 1-st and 2-nd order phase transitions as well as for the 
critical point p, Q] . 

In addition to the statistical fluctuations the complicated time evolution of A+A col- 
lisions generates dynamical fluctuations. The fluctuations in the initial energy deposited 
inelastically in the statistical system yield dynamical fluctuations of all macroscopic param- 
eters, like the total entropy or strangeness content. The observable consequences of the 
initial energy density fluctuations are sensitive to the equation of state, and can therefore 
be useful as signals for phase transitions [12]. Even when the data are obtained with a 
centrality trigger the number of nucleons participating in inelastic collisions still fluctuates 
considerably. In the language of statistical mechanics, these fluctuations in the participant 
nucleon number correspond to volume fluctuations. Secondary particle multiplicities scale 
linearly with the volume, hence, volume fluctuations translate directly to particle number 
fluctuations. 

The present work is a continuation of our recent study ^3] where we have analyzed the 
charged particle number fluctuations in Pb+Pb collisions at 158 AGeV within the UrQMD 
and HSD transport approaches. The net baryon number and electric charge event-by-event 
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fluctuations are studied in different rapidity regions of the projectile and target hemispheres. 

The paper is organized as follows. Section II presents the HSD results for the fluctua- 
tions of the number of nucleon participants while Sections III and IV give the net baryon 
number fluctuations and electric charge fluctuations, respectively. In Section V we discuss 
the fluctuations in the samples of most central collisions, Section VI shows a comparison of 
our calculations with experimental data from the NA49 Collaboration, whereas Section VII 
finally concludes the present study. 

II. FLUCTUATIONS OF THE NUMBER OF PARTICIPANTS 

In each A+A collision only a fraction of all 2 A nucleons interact. These are called par- 
ticipant nucleons and are denoted as Np^°-^ and Np^^ for the projectile and target nuclei, 
respectively. The nucleons, which do not interact, are called the projectile and target spec- 
tators, Ng^"-' = A — N^"-' and N^^^^ = A — N^p^^ . The fluctuations in high energy A+A 
collisions are dominated by a geometrical variation of the impact parameter. However, even 
for the fixed impact parameter the number of participants, Np = N^"-' + N'^^^ , fluctuates 
from event to event. This is due to the fluctuations of the initial states of the colliding 
nuclei and the probabilistic character of the interaction process. The fluctuations of Np 
form usually a large and uninteresting background. In order to minimize its contribution 
the NA49 Collaboration has selected samples of collisions with a fixed numbers of projectile 
participants. This selection is possible due to a measurement of N^^"-' in each individual 
collision by a calorimeter which covers the projectile fragmentation domain. However, even 
in the samples with N^"-' = const the number of target participants fluctuates considerably. 
Hence, an asymmetry between projectile and target participants is introduced, i.e. N^"-' is 
constant by constraint, whereas A'p '"^ fluctuates independently. 

In the following the variance, Var{n) = {n"^) — (n)^, and scaled variance, Un ^ 
Var{n) / (n), where n stands for a given random variable and (■ ■ ■ ) for event-by-event averag- 
ing, will be used to quantify fluctuations. In each sample with N^"-' = const the number of 
target participants fluctuates around its mean value, {N^p"^^) with the scaled variance u/p"^^ . 
From an output of the HSD minimum bias simulations of Pb+Pb collisions at 158 AGeV 
we form the samples of events with fixed values of N^°-' . Fig. ^ presents the HSD average 
value {N^p"^^) (left) and the scaled variances u'p^ (right) as functions of N^"\ One finds 
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{N*p^^) ^ N^"-'] the deviations are only seen at very small {N^°-' ^ 1) and very large 
{Np'°-^ ^ A) numbers of projectile participants. The fluctuations of N'^^^ are quite strong: 
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FIG. 1: The HSD simulations in Pb+Pb collisions at 158 AGeV for the average value {Np'^^) 
(left) and the scaled variances Wp**^ (right) as functions of N^"-' . 

The consequences of the asymmetry between projectile and target hemispheres depend 
on the A+A dynamics. According to Ref. different models of hadron production in 
relativistic A+A collisions can be divided into three limiting groups: transparency (T-), 
mixing (M-), and reflection (R-) models. The rapidity distributions resulting from the T-, 
M-, and R-models are sketched in Fig.|21taken from Ref. 3]. We note that there are models 
which assume the mixing of hadron production sources, however, the transparency of baryon 
flows, e.g. three-fluid hydrodynamical model Q]- R-models appear rather unrealistic and 
are included for completeness in our discussion. 



III. NET BARYON NUMBER FLUCTUATIONS 



We begin with a quantitative discussion by flrst considering the fluctuations of the net 
baryon number in different regions of the participant domain in collisions of two identical 
nuclei. These fluctuations are most closely related to the fluctuations of the number of 
participant nucleons because of baryon number conservation. 

The HSD results for ub in Pb+Pb at 158 AGeV are presented in Fig. IHl In each event 
we subtract the nucleon spectators when counting the number of baryons. The net baryon 
number in the full phase space, B = Ns—N-g, equals then to the total number of participants 
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FIG. 2: The sketch of the rapidity distributions of the baryon number or the particle production 
sources (horizontal rectangles) in nucleus-nucleus collisions resulting from the transparency, mixing 
and reflection models. The spectator nucleons are indicated by the vertical rectangles. In the 
collisions with a fixed number of projectile spectators only matter related to the target shows 
significant fluctuations (vertical arrows). See Ref. for more details. 

Np = Np^'^ + Np'"-' . At fixed Np'"-' the Np number fluctuates due to fluctuations of Np^^ . 
These fluctuations correspond to an average value, {Np'^^) ~ Np'"-' , and a scaled variance, 
uj^p^^ (see Fig. P). Thus, for the net baryon number fluctuations in the full phase space we 
flnd, 

_ Var{Np) {{Nr'Y) - {N^r _ 1 tar, 
- " ^{]^) - 2 • 

A factor 1/2 in the right hand side of Eq. (Q) appears because only half of the total number 

of participants fluctuates. 

Let us introduce and Up, where the superscripts p and t mark quantities measured 

in the projectile and target momentum hemispheres, respectively. Fig. El demonstrates that 
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FIG. 3: The HSD simulations for Pb+Pb collisions at 158 AGeV for fixed values of A^|,™^'. Left: 
The baryon number fluctuations in full acceptance, a;^, in projectile hemisphere, a;^ (lower curve), 
and in target hemisphere, o;^ (upper curve). The dashed line, 0.5 u/p^^ , demonstrates the validity 
of the relation Right: The scaled variances of the baryon number fluctuations in different 

rapidity intervals. 

cj^ > cu^, both in the whole projectile-target hemispheres and in the symmetric rapidity 
intervals. On the other hand one observes that u/^ ~ in most central collisions. This 
is because the fluctuations of the target participants become negligible in this case, i.e. 
^targ _^ g (^pjg^ [jj ^jght). As a conscqucncc the fluctuations of any observable in the 
symmetric rapidity intervals become identical in most central collisions. Note also that 
transparency-mixing effects are different at different rapidities. From Fig. ^ (right) it follows 
that cu^ in the target rapidity interval [—2, —1] is much larger than tu^ in the symmetric 
projectile rapidity interval [1,2]. This fact reveals the strong transparency effects. On the 
other hand, the behavior is different in symmetric rapidity intervals near the midrapidity. 
From Fig. ^ (right) one observes that uf^ in the target rapidity interval [—1, 0] is already 
much closer to in the symmetric projectile rapidity interval [0,1]. This gives a rough 
estimate of the width. Ay ^ 1, for the region in rapidity space where projectile and target 
nucleons communicate to each others. 

By assumption, the mixing of the projectile and target participants is absent in T- and 
R-models. Therefore, in T-models, the net baryon number in the projectile hemisphere 
equals to N^'^°^ and does not fluctuate, i.e. uj^{T) = 0, whereas the net baryon number in 
the target hemisphere equals to iV*"*"^ and fluctuates with to'^(T) = Up^^. These relations 
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are reversed in R-models. We introduce now a mixing of baryons between the projectile 
and target hemispheres. Let a be the probabihty for a (projectile) target participant to 
be detected in the (target) projectile hemisphere. We denote by n* and "nP the number of 
baryons which end uo in the target and projectile hemisphere, respectively, from the opposite 
hemisphere. Then the probabilities to detect 5* baryons in the target hemisphere, and B"^ 
baryons in the projectile hemisphere, can be written as. 



jqtarg j^proj 



ay p 



X - arr'-- S (B' - N'r - n' + W) , (2) 



nP\{Np'' - nP)\ 



J\ftarg n* = l nP = l 



'P 



X a-\l - afr^--' ^'^Z/ 6 (5- - - + , (3) 

n*!(A'p — n*)! 

where W{Np^^; Np'"-') is the probability distribution of Np^^ in a sample with fixed value 
of N^"-' . From Eqs. (j2l3|) with a straightforward calculation we find: 

Jp = {I - af + 2a(l - a) , = o? + 2a(l - a) . (4) 

A (complete) mixing of the projectile and target participants is assumed in M-models. Thus 
each participant nucleon with equal probability, a = 1/2, can be found either in the target 
or in projectile hemispheres. In M-models the fluctuations in both projectile and target 
hemispheres are identical. The limiting cases, a = and a = 1, of Eq. correspond to 
T- and R-models, respectively. In summary, the scaled variances of the net baryon number 
fluctuations in the projectile, Up, and target, Up, hemispheres are: 

= , coUT) = uj'r , (5) 

o;^(M) = Jp{M) = \+ \ , (6) 
u^piR) = Jp , Jp{R) = , (7) 

in the T- (0), M- © and R- ((Tj) models of the baryon number flow. The different models 
lead to significantly different predictions for cj^ and u^p. 
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FIG. 4: The a;^ {left) and lu^ (right) of the HSD simulations in comparison to T-, M- and R-models 
H5l7p , with uj^p^^ taken from Fig. ^ 

In Fig. |3]we show the predictions of T-, M- and R-models (jSJITj) with oo^p^^ from Fig. ^ 
for Pb+Pb collisions at 158A GeV. From Fig. 0] one concludes that the HSD results are 
close to the T-model estimates for baryon flow. However, the deviations from the results 
(0) are clearly seen: ujp>0 and ujp > u/p"^^ . One can not fit the HSD values of cj^ 
and uj'^p by Eq. (jl)). To make cu^ > one needs a > 0, but this induces uj^p < u/p"^^ , 
i.e. a mixing of baryons between the projectile and target hemispheres creates a non-zero 
baryon number fluctuations in the projectile hemisphere on the expense of fluctuations in 
the target hemisphere. Indeed, it follows from Eq. that u/p increases with a for all a, if 
^targ ^ a < (2 — cjp^^)^^, if u/p"^^ < 1. On the other hand, Up increases with a if 

a < (1 — uj*p^^){2 — cup^^)^-*^. This shows that an increase of uip with a is only possible for 
^targ ^ Thus for LUp^^ > 1 One finds an increase of Up with a and a decrease of ujp with 
a for all physical values of a from to 1. Therefore, we conclude that the HSD values of 
ujp (i.e. the fact that Up > ujp^^) can not be explained by Eq. Q with a > 0. 

The numbers of target and projectile participants are defined as N'^^^ = A — N^^^^ and 
jy-proj ^ A—Ng~°-^. The actual event-by-event numbers of baryons in the target and projectile 
hemispheres, Np and A'^^, may differ from Np^^ and N^"-' . This is because a transfer of 
baryons between the projectile and target hemispheres arises from the production of baryon- 
antibaryon pairs. The partners of each newly created 66-pair can be detected with non-zero 
probability in different hemispheres. We introduce 6* = Np — N^p^^ and the number of 
antibaryons in the target hemisphere, 6*. Similarly, V = Np — Np''"-' , while V is the number 
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of antibaryons in the projectile hemisphere. One finds: 



+ 



Var{N',^'3 + 6* - b ) 
W) 



targ 



Var{b') + Var{b') + 2 A(A^^'^"^ 6*) - 2 A{N'p^''' , b' 



2 A(bK b 



Var{N'';°' + y - b") 
W) 



where 



Var(lF) + Far(?) - 2 AiW , b 



(8) 

(9) 



A(iVi, N2) = (iVi . N2) - (iVi) ■ {N2) . 



(10) 



As N^°-' = const in the sample, it follows that 
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0, A{NP;"\ ¥>) = 0, A(iV^™^ fo^) = 



these terms are absent in the r.h.s. of Eq. 0. Different terms of Eq. (jH)) and Eq. 
found from the HSD simulations are presented in Fig |3| One observes that terms of 
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FIG. 5: Different terms of Eq. Q, left, and Eq. 0, right, are presented as a function of N^°-' . 

Eq. ()8I9|) expressing the fluctuations of antibaryons, Var(lj')/N^°-^ , and the correlation 
terms, 2A(A^p ^ F)/A^^™^ and - 2A(6*, b )/N'p;°\ with antibaryons included, are small. 
Therefore, one flnds, cj^ = Var{lP)/N^°-' . In the target hemisphere, the u/p^^ gives the 
main contribution to uj*q in Eq. (jH}. The term Var{b^) /N^°-' also contributes to tu^, sim- 
ilarly to that, Var{b^)/N^°-' , in the projectile hemisphere. However, the main additional 
term to tu^ is 2A{N'^^^, b^)/Np'^°-', which is due to (positive) correlations between N^^^ 
and 6*. This implies that in events with large N^p^^ (i.e. N^p^^ > {N^p^^) = N^°-^) some 
additional baryons move from the projectile to the target hemisphere, and when N'^'^^ is 
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FIG. 6: The HSD results for Pb+Pb collisions at 158 AGeV for the rapidity distributions of baryon 
numbers in nonsymmetric samples with N^"-' = 50,N^p^^ = 78 (left), and Np °-^ = 50, A'jf^ = 20 
(right). 

small (i.e. N^p^^ < {N^p'"^) = N^p"-') the baryons move in the reverse direction from the 
target to the projectile hemisphere as shown in Fig. IHl 

This HSD result looks rather unexpected. We remind that Eq. Q predicts for tu^ the 
opposite behavior: due to a simple mixing of baryons between the target and projectile 
hemispheres the initially large fluctuations, tUp*^^, are transformed into smaller ones, u^p. 
It seems that the origin of this effect is the following: For N^p^^ > N^"-' each projectile 
nucleon interacts, in average, more often than the target nucleon. The projectile participant 
loses then a larger part of its energy, and in the rapidity space its position becomes closer 
to Ucm. = than the position of target participants. This gives to projectile participants 
more chances to move due to further rescatterings from projectile to target hemisphere, in 
a comparison with target participants to move in the opposite direction. For N^p"^^ < N^"-' 
there is a reverse situation. This fact was not taken into account in Eqs. ()2f3j) where it has 
been assumed that the mixing probability a is the same for projectile and target participants, 
and independent of Np^^ . 

IV. NET ELECTRIC CHARGE FLUCTUATIONS 

The T-, M- and R-models give very different predictions for ujp and Up for the samples 
of events with fixed values of N^p"^ . Additional interesting correlations between the -B* and 
numbers, as those seen in the HSD simulations, can be expected. Unfortunately, they 
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may be difficult to test experimentally as an identification of protons and a measurement 
of neutrons in a large acceptance in a single event is difficult. Measurements of the charged 
particle multiplicity in a large acceptance can be performed with the existing detectors. 
In this section we consider the HSD results for the net electric charge, Q, fluctuations. 
As Q = OAB in the initial heavy nuclei one can naively expect that Q fluctuations are 
quite similar to B fluctuations. We stress, however, a principal difference between Q and 
B in relativistic A+A collisions. Fig. [7| demonstrates the rapidity distributions of the net 
baryon number, B = Nb — (left), and total number of baryons, Nb + (right), for 
different centralities in Pb+Pb collisions at 158 AGeV. One observes that both quantities 
are very close to each other; the y-dependence and absolute values are very close for B and 
Nb — N-g distributions. This is, of course, because the number of antibaryons is rather small, 
Nb < Nb- 
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FIG. 7: The HSD rapidity distributions in Pb+Pb collisions at 158 AGeV for the net baryon 
number, B = Nb — N^ i^sft), and total number of baryons, Nb + -/V-g- (right), at different N^"-' 
and in the minimum bias (m.b.) sample. 

Fig. IHl shows the same as Fig. [7|but for the electric charge Q = N+ — N_ (left), and 
total number of charged particles, Nch = N+ + N^ (right). The y-dependence of dQ/dy 
and dNch/dy is quite different. Besides, the absolute values of Nch are about 10 times larger 
than those of Q. This implies that Q <^ A^_|_ ~ A^_. 

In the previous section we have used the scaled variance to quantify the measure of 
the net baryon fluctuations. It appears to be a useful variable as ujb is straightforwardly 
connected to u^p^^ and due to the relatively small number of antibaryons. Fig. |H1 tells that 
ujq is a bad measure of the electric charge fluctuations in high energy A+A collisions. One 
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FIG. 8: The same as in Fig. [7|but for the electric charge Q = — i^eft), and total number 
of charged particles, Nch = + (right). 



observes that ujq = Var{Q)/ (Q) is much larger than 1 simply due to the small value of (Q) 
in a comparison with Nj^ and A^_. If the A+A collision energy increases, it follows, (Q) 0, 
and thus uq ^ oo. The same will happen with ub, too, at much larger energies. A useful 
measure of the net electric charge fluctuations is the quantity (see, e.g., [l^): 

A value of Xq can be easily calculated for the Boltzmann ideal gas in the grand canonical 
ensemble. In this case the number of negative and positive particles fluctuates according 
to the Poisson distribution (i.e. = = 1), and the correlation between A^_|_ and N_ 
are absent (i.e. (N^N^) = {N^){N^)), so that Xq = 1. On the other hand, the canonical 
ensemble formulation (i.e. when Q = const fixed exactly for all microscopic states of the 
system) leads to Xq = 0. Fig. IHl shows the results of the HSD simulations for the full 
acceptance, for the projectile and target hemispheres (left), and also for symmetric rapidity 
intervals in the c.m.s. (right). 

The Q fluctuation in the full acceptance is due to N'^^^ fluctuations. As Q = OAB in 
colliding (heavy) nuclei, one may expect Var{Q) = 0.16 Var{B). In addition, (Nch) = 4:{Np) 
at 158 AGeV, so that one estimates Xq = 0.04 ub for the fluctuations in the full phase space. 
The actual values of Xq presented in Fig. El (left) are about 3 times larger. This is because 
of Q fluctuations due to different event-by-event values of proton and neutron participants 
even in a sample with fixed values of N^"-' and N^p"^^ . 

From Fig. El (right) one sees only a tiny difference between the Xq values in the symmetric 
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FIG. 9: Left: The HSD simulations in Pb+Pb collisions at 158 AGeV for Xq at different values 
of N^"-' in the full acceptance (lower curve), for the projectile (middle curve) and target (upper 
curve) hemispheres. Right: The same, but for symmetric rapidity intervals in the c.m.s. 

rapidity intervals in the projectile and target hemispheres, and slightly stronger effects for 
the whole projectile and target hemispheres (Fig. IHl right). In fact, the fluctuations of 
and are very different in the projectile and target hemispheres, and the scaled variances 
uj*j^ and uj*_ have a very strong A^^j?,^ ""'-dependence. This is shown in Fig. obtained in our 
previous study jl7|. 
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FIG. 10: The HSD results for the scaled variances of negatively {left) and positively [right) charged 
hadrons in Pb+Pb collisions at 158 AGeV for the projectile (lower curves) and target (upper curves) 
hemispheres. 
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X, 



The Xq can be presented in two equivalent forms 

^ A(iV+,iV_) 



{N+) ^ (iV_) 



(iVc/^) 



2 cu. 



{Nch) 



(12) 



Eq. (fT^ is valid for any region of the phase space: full phase space, projectile or target 
hemisphere, etc. As seen from Fig. both uj\ and uj[_ are large and strongly N^"-'- 
dependent. This is not seen in Xq because of strong correlations between A^^ and N^_, i.e. 
the term 2 A{N^, N_)/ {Net) compensates and u!_ terms in Eq. ()12|). This is also seen 
from Fig. ^2 A cancellation of strong A^^'^°-'-dependence in the target hemisphere takes 




FIG. 11: The HSD results for the scaled variances of all charged hadrons, Wc/j, in Pb+Pb collisions 
at 158 AGeV for the projectile (lower curve) and target (upper curve) hemispheres. 

place between the sum of o;^ and tui terms of Eq. (fT^ . and the a;*;j-term. 

Fig.lT^shows a comparison of the HSD results for Xq with NA49 data in Pb+Pb collisions 
at 158 AGeV for the forward rapidity interval 1.1 < y < 2.6 inside the projectile hemisphere 
with additional py-filter imposed. As an illustration, the HSD results in the symmetric 
backward rapidity interval —2.6 < y < —1.1 (target hemisphere) are also included. One 
observes no difference between the Xq results for the NA49 acceptance in the projectile 
and target hemispheres. The HSD values for uj+, uj_, and Uch are rather different in the 
projectile and target hemispheres for the NA49 acceptance (see Figs. ITUl and ITTj) . This is 
not seen in Fig.^Jfor Xq. As explained above a cancellation between uj+, uj- and ujch terms 
take place in Eq. ()12j). In fact, NA49 did not perform the Xq measurements. The Xg-data 
(solid dots) presented in Fig. [121 are obtained from Eq. (|T2|l using the NA49 data for 
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FIG. 12: The HSD results for Xq for Pb+Pb collisions at 158 AGeV for the forward rapidity 
interval 1.1 < y < 2.6 inside the projectile hemisphere. The solid dots are the estimates obtained 



20| (the errorbars are not indicated here). For 



from Eq. ()12() using the NA49 experimental data 
illustration, the HSD results in the symmetric backward rapidity interval —2.6 < y < —1.1 (target 
hemisphere) are also presented. 



u!-, and Uch as well as (A^+), {N^), and (Nch) [20!]. Such a procedure leads, however, to 
very large errors for Xq (which are not indicated in Fig. IT^ which excludes any conclusion 
about the (dis) agreement of HSD results with NA49 data. 



V. FLUCTUATIONS IN MOST CENTRAL COLLISIONS 

In this section we consider the baryon number and electric charge fluctuations in the 
symmetric rapidity interval [—y,y] in the c.m.s. for the most central Pb+Pb events. We 
chose the sample of most central events by restricting the impact parameter to 6 < 2 fm. 
It gives about 2% most central Pb+Pb collisions from the whole minimum bias sample. 
Fig. ^1 shows the HSD results for electric charge fluctuations in 2% most central Pb+Pb 
collisions for the symmetric rapidity interval AY = [—y, y] in the c.m.s. as the function 
of Ay = AY/2. For AY one finds Xq — > 1. This can be understood as follows: For 
Ay the fluctuations of negatively, positively and all charged particles behave as for 
the Poisson distribution: uj+ = = Uch = 1- Then from Eq. (fT^ it follows that Xq = 1, 
too. From Fig. El (right) one observes that uj+, uj^, and Uch all increase with increasing 
interval AY. However, Xq decreases with AY and - because of global Q conservation - it 



15 



Ql . , . , . 1 

12 3 

Ay 



FIG. 13: The HSD results for electric charge fluctuations in 2% most central Pb+Pb collisions at 
158 AGeV in the symmetric rapidity interval , Ay = [—y,y] as a function of Ay = AY/2 in the 
c.m.s. A left panel shows the behavior of Xq, and a right one demonstrates separately uj-^-, uj^, 
and ujch- 

goes approximately to zero when all final particles are accepted. 

In Fig. HH (left) the HSD results for the scaled variances are presented in full acceptance 
as functions of Np''"-' . FigOl (right) demonstrates the probability distribution of events with 
6 < 2 fm over N^°-' . One observes that even in the 2% centrality sample the values of A'^^'^"-' 




FIG. 14: The HSD results in Pb+Pb collisions at 158 AGeV. Left: The scaled variances a;+, 
and LiJc/i in the full acceptance. Right: The distributions of events over N^"-' in most central 
collisions with 6 < 2 fm. 

are noticeably smaller than the maximum value, A = 208. As seen from Fig. El (left) the 
HSD values of 6^;+, uj^, and uJch become then essentially larger than 1 in agreement with 
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those presented in Fig. [121 

Fig. El shows the net baryon number fluctuations in the symmetric rapidity interval 
[— y,?/] in the c.m.s. as the function of AY. As a measure of the net baryon number 



— ■ — antibaryons 
— — baryons 
—•—total 




FIG. 15: The HSD results for net baryon number fluctuations in 2% most central Pb+Pb collisions 
at 158 AGeV in the symmetric rapidity interval Ay = [— y, y] as a function of Ay = Ay/2 in the 
c.m.s. The left panel shows the behavior of Xb, and a right panel presents separately ojnb^ '^^b' 
and WAT^+AT . 



fluctuations we have used the quantity, 

Xb = 



Var{B) 



(13) 



{Nb + N^) 

As for the electric charge, one flnds that Xb ^ 1 at /XY (this is because all ojnbi ^n^i 
and ujNg^N-g go to 1 in this limit (see Fig. UHl left), and Xb at upper limit of XY 
because of global baryon number conservation. 
Writing the variance Var{B) in the form. 



Var{B) = 2 Var{NB) + 2 Var{N^) - Var{NB + N-^) 



we flnd 



X. 



2 ujnp 



(Nb) 



+ 2 UN-, 



{Nb) 



(14) 



(15) 



{Nb + N^^ {Nb + N^) 

The behavior of the different terms in Eq. (fT3|) is the following: As seen from Fig. ^1 right, 
cuat— = 1 for all values of AY. This is because N^ <^ Nb, and baryon number conservation 
does not affect the fluctuations of antibaryons. Due to the small number of antibaryons in 
comparison to baryons, one also observes ub = ^^Nb — ^Nb+n— 
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VI. ELECTRIC CHARGE FLUCTUATIONS IN CENTRAL PB+PB COLLISIONS 
AT 20, 30, 40, 80 AND 160 A GEV 



In this section we present the HSD results for the event-by-event electric charge fluctu- 
ations as measured by the NA49 Collaboration in central Pb+Pb collisions at 20, 30, 40, 
80 and 160 A GeV 22j. The interest in this observable (as a signal of deconfinement) is 
related to the predicted in Refs. [3, suppression of event-by-event fluctuations of the 
electric charge in a quark-gluon plasma relative to a hadron gas. However, these predic- 
tions were based on the assumption that the initial electric charge fluctuations survive the 
hadronization phase. 

The first experimental measurement of charge fluctuations in central heavy-ion collisions 



22j at SPS showed a quite mod- 



by PHENIX |25] and STAR |26| at RHIC and by the NA49 
erate suppression of the electric charge fluctuations. This observation has been attributed 
to the fact that the initial fluctuations are distorted by the hadronization. In particular, the 
observed fluctuations might be related to the final resonance decays. 

In this respect it is important to compare the experimental data with the results of mi- 



croscopic transport models such as HSD where the resonance decays are included by default. 
In order to quantify the event- 
quantity $ defined as 



In order to quantify the event-by-event electric charge fluctuations we have calculated the 



^ = ^/^-V^^ (16) 



where 

TV 



q-q, Z = 5^(g,-g). (17) 



i=l 



Here q denotes a single particle variable, i.e. electric charge g; N is the number of particles of 
the event within the acceptance, and over-line and (...) denote averaging over a single particle 
inclusive distribution and over events, respectively. By construction, $ of the system, which 
is an independent sum of identical sources of particles, is equal to the $ for a single source 
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28|. 



In order to remove the sensitivity of the final signal to the trivial global charge conserva- 
tion (GCC) the measure A$q is defined as the difference: 

= $,,Gcc • (18) 
18 



Here the value of $g is given by j29|, |30| : 



g,GCC 



(19) 



where 



p = -f^ (20) 

V^ch/tot 

with {Nch) and {Nch)tot being the mean charged multiphcity in the detector acceptance and 
in full phase space (excluding spectator nucleons), respectively. 

By construction, the value of A$q is zero if the particles are correlated by global charge 
conservation only. It is negative in case of an additional correlation between positively 
and negatively charged particles, and it is positive if the positive and negative particles are 
anti-correlated js^. 



Energy 

[A GeV] NA49 HSD 
20 ■ — Q— 
30 
40 



<N >/<N 
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FIG. 16: The dependence of the $q (l.h.s.) and A<I>g (r.h.s.) on the fraction of accepted particles 
for central Pb+Pb collisions at 20-158 ^GeV. The NA49 data are shown as full symbols, 
whereas the open symbols (connected by lines) stay for the HSD results. The dashed line shows 
the dependence expected for the case if the only source of particle correlations is the global charge 
conservation ^q^ccCi Eg. [TUl 

Figure ITHl shows the HSD results for the dependence of $q (l.h.s.) and A$g (r.h.s.) on 
the fraction of accepted particles (Nch) and {Nch)tot (calculated for ten different rapidity 
intervals increasing in size from Ay = 0.3 to Ay = 3 in equal steps) for central Pb+Pb 
collisions at 20, 30, 40, 80 and 158 A GeV. The NA49 data ^ are shown as fuU symbols. 
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whereas the open symbols (connected by hnes) reflect the HSD results. The dashed line 
shows the dependence expected for the case if the only source of particle correlations is the 
global charge conservation $g,GCC (Eq. (HSI))- 

The data as well as the HSD results for $q (Fig. E[ l.h.s.) are in a good agreement and 
show a monotonic decrease with increasing fraction of accepted particles. After substraction 
the contribution by global charge conservation (the dashed line in Fig. IT^ . the values of 
A*, va.y between and -0.05 which „e .igniflcantly la.ge. than the vainer expected tor 
QGP fluctuations (-0.5 < A<l>q < -0.15 |30|). 




FIG. 17: The energy dependence of A<I>g measured in central Pb+Pb collisions for a narrow rapidity 
interval Ay = 1.2 (l.h.s.) and a broad rapidity interval Ay = 3 (r.h.s.). The NA49 data are 
shown as full symbols, whereas the the open symbols (connected by lines) reflect the HSD results. 

Figure El presents the energy dependences of for two selected rapidity intervals - the 
intermediate rapidity interval Ay = 1.2 (l.h.s.) and for the largest rapidity interval Ay = 3 
(r.h.s.). The both, data and HSD results, show the a weak decrease of A$g with increasing 
energy. 

The fact that the HSD model, that includes no explicit phase transition, describes the 
experimental data can be considered an independent proof that the event-by-event charge 
fluctuations are driven by the hadronization phase and dominantly by the resonance decays 
(which are naturally included in HSD) and no longer sensitive to the initial phase fluctuations 
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from a QGP. 



VII. SUMMARY AND CONCLUSIONS 

The goal of this study was to investigate the sensitivity of event-by-event fluctuations 
of baryon number and electric charge to the early stage dynamics of hot and dence nu- 
clear matter created in heavy-ion collisions at SPS energies and the infiuance of the futher 
hadronization and rescattering phase. For that perpose we have explored the microscopic 
HSD transport model which allows also to investigate (on event-by-event basis) the influence 
of the experimental acceptance and the set-up on the final observables. 

It has been found that the fluctuations in the number of target participants strongly 
influences the baryon number and charged multiplicity fluctuations. The consequences of 
this fact depend crucially on the dynamics of the initial flows of the conserved charges and 
inelastic energy. 

For a better quantitative understanding of the microscopic transport model (HSD) results 
we have considered 3 limiting groups of models for nucleus-nucleus collisions: transparency, 
mixing and reflection. These "pedagogical" considerations indicate that the HSD model 
(as well as UrQMD, cf. Ref. |l3]) shows only a small mixing on initial baryon flow and is 
closer to the T-model. This supports the findings from Ref. (2] about the influence of the 
partonic degrees of freedom on the initial phase dynamics which might increase the mixing 
by additional strong parton-parton interactions. Thus, the measurement of the net baryon 
number fluctuations helps to quantify the mixing of initial baryon flow. 

The first microscopic event-by-event calculations of the charge fluctuations A$g within 
the HSD model show a good agreement with the NA49 data at SPS energies. Thus, this 
observable is dominated by the final stage danymics, i.e. the hadronization phase and the 
resonance decays, and rather insensitive to the initial QGP dynamics. 
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